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We theoretically analyze the dynamical evolution of photonic quantum walks on Mobius strips and 
other exotic structures in 3D integrated photonics. Our flexible design allows discrete observations 
of continuous time quantum walks of photons in a variety of waveguide arrays. Furthermore, our 
design allows one to inject photons during the evolution, allowing the possibility of interacting 
with the photons as they are 'walking'. We find that non-trivial array topologies introduce novel 
time-dependent symmetries of the two-photon correlations. These properties allow a large degree 
of control for quantum state engineering of multimode entangled states in these devices. 



I. INTRODUCTION 

With current advances in quantum technologies, quan- 
tum walks are beginning to be experimentally explored 
in a variety of systems. The experimental study of 
quantum walks is of practical importance in quantum 
simulation [T] and solving hard computational problems 
including the boson sampling problem [2] and graph 
searches [3]. Recently, quantum walks have been ob- 
served in waveguide arrays 0HS], fiber loops [7J [8] and 
NMR [9 . Furthermore, quantum walks have been pro- 
posed in several other systems including ion traps [10] 
and Rydberg atoms [TTj. 

In particular, waveguide arrays in integrated photonics 
are a promising experimental implementation of quan- 
tum walks due to their flexible geometry and low losses. 
Recent experiments in integrated photonics have demon- 
strated a continuous time quantum walk (CTQW) in a 
ID waveguide array [1] and a discrete time quantum 
walk in a quasi- ID array of directional couplers [5]- A 
limitation of these approaches is that the length of the 
waveguide array is fixed and therefore it is only possible 
to observe the walk at a single time instant. Although 
it is possible, in principle, to reduce the array length by 
cut-back techniques, this 'ad hoc' process is highly time 
consuming and achieving consistent coupling at each cut- 
back stage to allow meaningful comparison of photon cor- 
relations would be very tough. 

Here, we describe an experimental proposal for a 3D 
waveguide array implementation of a CTQW, which can 
be observed at multiple discrete times (or discrete optical 
path lengths depending on the reader's preferred mental 
picture) . This generalizes previous waveguide approaches 
which could only observe CTQWs at a single time in- 
stant. The approach allows a variety of exotic structures 
with varying boundary conditions to be studied. An- 
other novel aspect of our proposal is that photons can be 
added to the CTQW during its evolution. We show that 
by varying the waveguide boundary conditions and the 
time delay between the injection of photons, it is possible 
to engineer and analyze a range of multimode quantum 
states with time-dependent photon statistics. 



This paper is structured as follows, in section[TT]we dis- 
cuss our proposal for a general waveguide array. We then 
calculate the experimentally relevant two-photon corre- 
lation function for simultaneous and delayed two-photon 
input states in section [TTTl We next analyze the dynamics 



of three particular examples of our proposal in section IV 



a cylindrical array, a Mobius strip waveguide array and a 
twisted circular waveguide array. In section [V] we discuss 
practical aspects of an experimental implementation. We 
conclude in section IVll 



II. GENERAL DESCRIPTION OF THE DEVICE 

We begin with the dynamics of a general implemen- 
tation of our waveguide array device. We consider the 
concep tual d evice depicted in Fig. |l(a)| with three exam- 
ples in l(b)| -(d). Initially, photons are injected into a set 
of directional couplers which either couple into a waveg- 
uide array or continue into a set of output modes to be 
detected. We then post-select events where all photons 
enter the waveguide array by monitoring the arrival times 
of the output photons. We describe the meaning of the 
generic operators B, H and P shortly. 

Three particular examples of this general waveguide 
device are presented in Fig. 1(b) -(d): the cylindrical ar- 
ray, the Mobius strip array and the twisted circular ar- 
ray. Each device has a radius much larger than the inter- 
waveguide distance and an array of directional couplers 
to couple the light in and out of the device. At the cou- 
pling region the waveguides fan-out so that each direc- 
tional coupler only couples to a single waveguide. This 
region can also be used to correct any path length dif- 
ferences between the waveguides that is acquired during 
the loop. 

The dynamics of the photons in the general waveguide 
array are determined by the Hamiltonian, 
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where, G n ,m is the rate of coupling between the n-th 
and m-th waveguide. As the coupling matrix is real and 
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FIG. 1. a) General schematic of experimental proposal. Photons are injected into directional couplers (denoted by B), which 
couple into a waveguide array and evolve under a local Hamiltonian H for time r. Due to the geometry of the array, the 
waveguides then undergo a permutation denoted by P. After the permutation the photons enter the directional couplers again 
and either exit the device to be detected or remain in the waveguide array for another loop. b)-d) Three specific implementation 
of the proposal: b) cylindrical array c) Mobius strip waveguide array and d) twisted circular waveguide array. For clarity the 
directional coupler array on the twisted circular array has been omitted. 



Hermitian, G n . m = G m , ra , it is possible to decompose it 
as G = VNV\ where Aij = 5ij\j is a diagonal matrix of 
eigenvalues and the j-th column of the unitary matrix, V, 
is the j-th eigenvector of G. We denote the n-th element 
of the m-th eigenvector of G as v n . m . Moving to the 
Heisenberg picture and using the eigendecomposition we 
can find the evolution of the operators [H] , 



(*) = E£Li Uj, k (t)& k (o), 



where 
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These relations are sufficient to find the evolution of an 
arbitrary state in the waveguide array. 

The input photons evolve in the waveguide array under 
the Hamiltonian for time t = L/v g , where r is the 
time it takes for a single traversal of the waveguide array 
loop, L is the path length of the loop and v g is the group 
velocity of photons in the array. (The structure is as- 
sumed to be engineered such that the group velocity and 
path length of each loop are the same) . After the system 



evolves for r, due to the topology of the array the loop 
waveguide modes as labeled by the input/output cou- 
plers have encountered an effective permutation which 
we make explicit. For example, in the Mobius strip array 



depicted in Fig. 1(c) the permutation relabels the modes, 
aj — > aN+i-j, whereas for the twisted circular array in 

%j+ c f° r some integer 



1(d) the permutation is a 7 - — > a 



Fig 

c. The cylindrical array in Fig. 1(b) induces just the iden- 
tity permutation. In general we denote the permutation 
by the operator P which has the action, 



a, — > a 



p(i)' 



(4) 



where p(j) is a bijective function of the mode num- 
ber. We denote n compositions of the permutation 
function as, p„(j). For example, p 3 (j) = p(p(p(j))), 
and similarly for the inverse function, p^O) = 

p-^p-^O'))) ■ 

After the permutation, the photons encounter the di- 
rectional couplers (the waveguide analog of a beamsplit- 
ter). The couplers act on the waveguides in the array 
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(a,,) and the input /output waveguides (bj) as 

a] -> cos(^j)a] + i sin(0j)&t, (5a) 
S] ->■ i sin(^)a] + cos(0j)St. (5b) 

With some probability (determined by 6j) a number 
of photons may exit the device to be detected. Alter- 
natively, there is some probability that all photons will 
remain in the device. In this case, the photons will again 
evolve under H for time r, then undergo the permuta- 
tion P before interacting once again with the directional 
couplers. Depending on the reflectivities of the couplers, 
all photons may remain in the waveguide array for a sub- 
stantial period of time, t r. 

To summarize, photons enter the device via the input 
ports, then evolve under a local Hamiltonian; after this 
the waveguides undergo a permutation then interact with 
a set of directional couplers. Any photons that exit the 
device are measured and the photons remaining in the 
device repeat the process. The projected wavefunction 
after n traversals of the device is 

n 

\iP(nr)) = [] MBPe-i" T mO)), (6) 

k=l 

where B denotes the action of the N directional couplers 
and M denotes a projective measurement on the output 
bj modes. 

The interaction with the directional couplers after each 
loop allows us to observe the Hamiltonian evolution at 
discrete time periods. Therefore the device can be used 
to discretely observe CTQWs of photons in a range of 
waveguide arrays. In addition, the permutation of the 
waveguides each step causes the device to exhibit inter- 
esting time-dependent symmetries of the photon statis- 
tics. By varying the initial state and selecting particular 
evolution times the device can be used to prepare a range 
of multimode entangled output states with desired sym- 
metry properties. 

Finally we note that there is a class of states that are 
simultaneous eigenstates of both P and H, To find these 
states we first write the Hamiltonian as a sum of normal 
modes, H = ^J]„ =1 A„4c„, where, dj = £) fc= i Uj.kOfc- 
Eigenstates of H which are invariant under the permuta- 
tion have the property P£j — £j , from which we find the 
condition Ujp-i(fc) — v j,k- In particular, two-photon si- 
multaneous eigenstates have the interesting property that 
the observed two-photon correlation function is identical 
for all time steps, despite the permutation of the waveg- 
uides and the Hamiltonian evolution. 



III. TWO-PHOTON INPUT 

We now illustrate the evolution of photons in the de- 
vice by giving two experimentally relevant examples. We 
first consider injecting two-photons into the device simul- 
taneously for a CTQW and observing the time-dependent 



photon statistics. We also consider delaying the input of 
the second photon, which allows the experimentalist to 
interact with the walk as it evolves. 



A. Simultaneous Two-Photon Input 

Here we consider injecting a single photon into in- 
put waveguides j and k and post-selecting the events 
where both photons enter the waveguide array. The post- 
selection is achieved by rejecting events where a photon is 
detected in an output waveguide for t < L/v g . If we fur- 
ther assume each directional coupler has the same angle 
6, the post-selected state occurs with probability sin 6. 
The post-selected initial state of the waveguide array is 
therefore 

Mo)> = a}o£|o>, (?) 

where the multimode vacuum is denoted |0) = 

io>no>r- 

The state then evolves under the Hamiltonian H for 
time r. After this evolution and a permutation of the 
modes, the state becomes 

mr))=Pat(-r)a{(-r)\0) 

N 

= E ^,m(r)^, g (r)at (TO) at (g) |o), 
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where the evolution of the modes, fit(— t), was found us- 
ing ([2]) and the relation Q was used to permute the 
modes. 

This state then interacts with N directional couplers 
according to |5]), 

N 

W T )) = E U 3,p- 1 {r){r)U k , p -x(s){T) 

x ^al cos 9 + ibl sin Q^j 

x faJcos0-Mbisin6>) |0), (8) 

where we have relabeled the summation indices to make 
explicit the permutations' effect on the coefficient of the 
mode operators. From this expression we find that at 
time r, the probability that both photons exit the device 
is sin 4 6 and the probability that both photons remain in 
the device is cos 4 9. 

If both photons exit the device we can detect signa- 
tures of quantum interference using the two-photon cor- 
relation function, T{' k s (t) HH3Q3]. This is defined as 
the probability that a single photon is detected in out- 
put waveguides r and s at time t, given that initially a 
single photon was injected into waveguides j and k, 

r#(r) - n ^-(V'(r)|St(0)St(0)S.(0)6 r (0)|^(r)).(9) 

1 + <V,S 

Using the wavefunction ([8|, we find the two-photon cor- 
relation function after one traversal of the device is 
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sin 4 (0) 
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(10) 



We see that the two-photon correlation function is strongly dependent on the inverse permutation functions p _1 (7') 
and p _1 (s). Therefore the permutation is readily observable in the photon statistics. Furthermore, by comparing the 
two-photon correlation function to the classical correlation function for distinguishable particles [13] , 



P nh T ) = YTJ^ (l f7 ^p- 1 (r)(' r ) f/ 'fe,p- i w( r )| 2 + \ U 3,V-H>)( T ) U k,p- 1 (r)( T )\ 2 ) 



(11) 



we see that quantum interference arises due to the fact that the two terms in ( 10 ) are added prior to taking the 
modulus. 

Alternatively, if both photons remain in the waveguide array, then we can determine their state at time t = 2t 
by repeating the above steps of Hamiltonian evolution, permutation and the interaction with directional couplers. 
Continuing, for n time steps, t — nr, we find the two-photon correlation function, 



r#(nr) = 



cos 4(n-l)(g) sm 4(g) 
1 + S r . s 



U J, P 7 l 1 (r)( nT ) U k,p- 1 (s)( nT ) + U j,p- 1 (s)( nT ) U k^- 1 (r)( nT ) 



(12) 



We see that the two-photon correlation function at t = 
nr depends on n applications of the inverse permuta- 
tion function. This property leads to very interesting 
time-dependent photon statistics. For example for the 
Mobius strip considered in section [IVB the permutation 
function has the property, P2 1 (s) = s, and we find dif- 
ferent statistics at odd and even time periods. By gating 
the output modes to only allow output at a particular 
time step, it is possible to prepare a range of multimode 
entangled states. 



Another inte rest ing property of the two-photon corre- 
lation function { 12 1 is the pre-factor, cos 4 ( n_1 )((9) sin 4 (6'). 
This pre-factor reduces the probability of a coincidence 
detection at later times, due to the interaction with the 
directional couplers at previous time steps. As choos- 
ing a weaker interaction with the directional couplers in- 
creases the probability that both photons remain in the 
device for longer [T3], it is possible to optimize the cou- 
pling angle, 9, to maximize the two-photon correlation 
function at time t — nr. For a desired time t = tit, 
we find that cos 4 ( n_1 )(0) siri 4 (0) attains a maximum at 
cos(#) = \/n — 1/y/ri. Therefore, the two-photon cor- 
relation function r^'(nr) is maximized by choosing the 
optimal coupling angle, 9 op t — cos~ 1 (y/ri^T/y/n). By 



building a variety of devices with differing coupling an- 



gles, it is possible to observe CTQWs over a range of time 
periods. 



B. Two-Photon Input with Delay 

Our proposal also allows the possibility of injecting 
photons at different times. It is therefore possible to 
effectively interact with the quantum walk as it evolves. 
This novel feature is not possible in existing quantum 
walk experiments in integrated photonics. 

Here, we consider injecting a photon into waveguide j 
at t = and another into waveguide k at t = td = ridT, 
where rid is an integer. We then monitor the output 
waveguides and post select events where both photons 
remain in the array for t > td- The wavefunction of the 
photons in the array immediately after the injection of 
the second photon is, 

m > t d )) = a{(-(t - t d ))&l ndU) (-t)\0), (13) 

where the first photon has undergone rid permutations. 
We see that the second photon interacts with the per- 
muted single photon quantum walk after it has evolved 
for t d . 

Repeating the analysis of the previous section we find 
the two-photon correlation function at t — nr + ridT is 



rii(nT + n d T)=Af 



U A „-i 



3P 



(r)(( n + n d)T)U kp -l (s) {nT)+U ]p -l^^ 



(14) 



where N = cos 4 ^" 1 )^) sin 4 (0)/(l + 8 r>a ). The two- 
photon correlation function exhibits clear interference be- 
tween the two walkers at different times. This property 



I 

can create a range of entangled states, which can be found 
in a single experiment by varying the input ports and de- 
lay time td- 
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IV. EXAMPLES: TWO-PHOTON INPUT 

In this section we discuss three particular waveguide 
implementations of our proposal: the cylindrical array, 
the Mobius strip array and the twisted circular array. 
For each implementation we diagonalize the local Hamil- 
tonian, describe the geometry dependent permutation 
function and discuss the dynamics of a two-photon in- 
put state. 



A. Cylindrical Array 

The simplest implementation of our proposal is the 
cylindrical array presented in Fig. 1(b) Here, an equi- 
spaced one-dimensional (ID) 



coupling matrix (16). The eigenvalues and eigenvectors 



waveguide array is curved 
around the surface of a cylinder with a large radius. We 
require the radius to be the largest length scale in the sys- 
tem, much larger than both the inter-waveguide spacing 
and the coherence length of the injected photons. Pho- 
tons are injected into the input waveguides, which couple 
into the array via a series of directional couplers. The 
photons then evolve on the effective ID waveguide array 
around the surface of the cylinder. After one traversal 
of the cylinder, the photons interact with the directional 
couplers where some photons may couple into the output 
modes to be detected. The remaining photons will then 
traverse the cylinder a second time and interact with the 
couplers. Depending on the coupling angle 6, this pro- 
cess may occur for a large number of traversals of the 
cylinder. The relevant permutation is just the identity, 
P(j) =3- 

The equispaced ID waveguide array has recently been 
the subject of intense theoretical jTHj and experimental 
work [4]. In particular, as the photon detectors were 
placed at the end of the ID array, the experiment was 
only able to measure the array at one particular interac- 
tion time 0]. In contrast the implementation presented 
here allows the ID array to be observed at several dis- 
crete times with the same device. This ability allows the 
experiment to observe time-dependent properties of the 
CTQW. 

We now consider the evolution of the photons in the 
cylindrical waveguide array. The effective local Hamilto- 
nian of the array in Fig. 1(b) is the ID waveguide array 
Hamiltonian, 



N N-1 



3 ^ 



(15) 



where all waveguides have the same effective frequency w 
and couple to their nearest neighbors at the rate g. This 
Hamiltonian is a special case of ([!]) with 



G jt k = u)5j t k + g{5 jik+ i + Sj + i t k). 



(16) 



The evolution of the photons in the array can be deter- 
mined from the Heisenberg picture mode operators ([2]), 
which depend on the eigenvalues and eigenvectors of the 



of this matrix are well known 1151 



Aj : = uj + 2g cos 



N + l 



_yK_ 
N + l 

jkir 



N + l 



(17a) 
(17b) 



By substituting these relations into (|6]), we can find 
the evolution of the wavefunction of an arbitrary ini- 
tial state in the device. Furthermore, using (17) it 



is straightforward to derive the two-photon correlation 
function r^(nr) (12). The experimentally observable 
two-photon correlation function is shown in Fig. 2(a)[ - 
2(d) for several time steps. By injecting pairs of pho- 



tons and detecting coincidences over many experimental 
runs, it is possible to build up several time steps of the 
two-photon correlation function. This is in contrast to 
current experiments, which can only measure the two- 
photon correlation function at a single time step [3]. 

We also consider the case of two-photon input with 
delay in the cylindrical waveguide array. The time evo- 
lution of the two-photon correlation functions for three 
different delays are shown in Fig. [3} The top row in Fig. [3] 
has simultaneous two-photon input, the second row has 
the second photon delayed by r, and the third row has 
the second photon delayed by 2r. We see that states 
with different input delays have quite different statistics. 
Therefore by varying the delay, a variety of multimodc 
entangled states can be prepared. 



B. Mobius Strip Array 

We now turn to the Mobius strip array depicted in 
Fig. 1(c) This device is intrinsically three dimensional, 
and outside the domain of conventional lithography. Sim- 
ilar to the cylindrical waveguide array, photons are in- 
jected in the input modes, which couple into the array 
via a series of directional couplers. The waveguide array 
then traces out a large loop, with a single twist, forming 
a Mobius strip. The twist in the loop implies that af- 
ter a 2tt rotation, the first waveguide becomes the N-th 
waveguide and vice-versa. This geometrically induced 
permutation is described by the permutation function 
p(j) — N + 1 — j. This permutation has the property 
P2CO = j, which implies that after a Att rotation, the 
j-th waveguide is mapped onto itself. 

Assuming the radius of the strip is large, we can ap- 
proximate the local Hamiltonian of the array as the ID 



waveguide array Hamiltonian (15). Therefore the eigen- 



values and eigenvectors for the Mobius strip array are 
the same as the previous section. This approximation is 
only valid where the path lengths of the waveguides are 
identical, which is only approximately true for a strip of 
large radius. Experimentally, the path lengths will need 
to be adjusted to make the path lengths identical, as will 
be discussed in section Ivl 
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FIG. 2. Two-photon correlation functions r^'J(nr) (Eq. ( 12 1) for cylindrical array (top row) a) n = 0, b) n = 1, 
d) n — 3 and for the Mobius strip array (bottom row) e) n = 0, f) n = 
time in units of gr _1 and for clarity we re-scale the vertical axis by cos 4 ^™" 
statistics at odd time steps in the Mobius strip array on the bottom row. 



j) n = 2, 

1, g) n = 2, h) n = 3. Here we set r = 1, measure 
1) (0)sin 4 (6'). Note the mirror flipping of the photon 



The Mobius strip device has a novel topology for 
CTQWs, whereby the walkers repeatedly evolve for a pe- 
riod of time then undergo a permutation. This topology 
is readily observable in the photon statistics, through the 
dependence of the correlation function ( 12 ) on the per- 
mutation function. The two-photon correlation for the 



Mobius strip is shown in Fig. 2(c) - 2(h)) where we see 
that the strip topology induces a mirror reflection of the 
statistics at odd time steps. Therefore, this device al- 
lows discrete time observation of a CTQW with time- 
dependent symmetries of the photon statistics. 



C. Twisted Circular Array 

We now consider the twisted circular array depicted 
This array is a 3D generalization of the 2D 



in Fig. 1(d) 



waveguide array recently experimentally studied [6] . The 
array is formed by bending a 2D circular waveguide array 
into a large loop in the shape of a torus. By rotating the 
waveguides in the 2D array during the loop, it is possible 
to introduce a permutation of the modes. For example 
for N — 3 waveguides, a 27r/3 rotation during the loop 
causes the permutation, a\ — > a-2, &2 — > 03, S3 — > a%. 

The local Hamiltonian of the array, including the cou- 
pling between all non-adjacent modes, is, 



N N 



(18) 



where due to the circular geometry, the subscript, n + a — 
1, is modulo N, and due to symmetry, gj > N_+% = 9N—j+2- 



Here the coupling matrix, G, is a circulant matrix defined 
by the vector (51,32, ■ ■ ■ ,9n)- For example for N = 5, 



G = 



(9x 92 33 34 35 \ 

35 31 92 33 34 

34 35 31 92 33 

33 34 35 3i 32 

V32 33 34 35 9iJ 



(9l 92 33 33 32 \ 

32 31 92 33 33 

33 32 31 92 33 
33 33 32 31 92 

V32 33 33 32 9iJ 



(19) 



As G is a circulant matrix, it can be diagonalized by the 
discrete Fourier transform (DFT) . Using the DFT we find 
the eigenvalues and eigenvectors of G are, 



(j-i)(/c-i) 



fc=i 
1 



N 



e N 



(20a) 
(20b) 



To our knowledge this is the first exact diagonalization of 
the circular waveguide array with non-nearest neighbor 
couplings. 

The twisting of the array as it completes a loop causes 
a permutation of the modes. The simplest permutation is 
a 2ir/N rotation of the 2D circular array as it completes 
as single loop, 



P(i) = 



j + l itj^N 
1 ifj = JV. 



(21) 



After N traversals of the loop, the modes return to their 
original order, PjvC?) — P(i)- The permutations from 
twisting the array leads to non-trivial photon statistics. 
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FIG. 3. Two-photon correlation functions r r ; 3 (w + n^r) (Eq. ( 14 () for cylindrical array with delayed two-photon input. Top 
row: rid — (simultaneous two-photon input) a) n = b) n = 1 c) n = 2 d) n — 3. Middle row: rid = 1 (input the second 
photon after one traversal) e) n — f) n — 1 g) n — 2 h) n = 3. Bottom row: rid = 2 (input the second photon after two 
traversals) i)n = 0j)n = lk)n = 2 1)n = 3 Here we set r = 1, measure time in units of g -1 and for clarity we re-scale the 
vertical axis by cos 4(n - 1) (6»)sin 4 (e). 



For example for N = 3, a 2tt/3 rotation causes the cor- 
relation matrix to permute as, 




(22) 



We see that all elements of the correlation matrix shift 
across and down one element, except those in the last 
row and last column which move to the first row and 
column. More complex permutations can be generated 
by multiple applications of this 2tt/N rotation of the 2D 
circular array each loop. 

The time evolution of the two-photon correlation ma- 
trix of a twisted circular waveguide array with twelve 
waveguides is shown in Fig.|4j The array has a 4 x 27r/12 
rotation each loop, which causes the permutation hj — ¥ 
a,j + 4 L . Hence after three traversals of the torus, the waveg- 



uides return to their original ordering. The effect of 
the permutation on the twisted circular array correla- 
tion function in Fig. [4] is more complex than those of the 
Mobius strip array in Fig. [2j Here, the elements of the 
correlation matrix shift across and down four elements 
each time step, except the last four rows and columns, 
which wrap around to to the first four rows and columns. 
By changing the number of 2tt/N rotations of the 2D cir- 
cular array each loop, a variety of multimode entangled 
states can be prepared. 



V. EXPERIMENTAL ISSUES 

We now consider the conditions required for experi- 
mental demonstration of our ideas. We consider two pos- 
sible implementations: waveguide structures written in 
glass blocks (typically boro-silicates) by the femtosecond 
laser direct-write (FLDW) technique [16] , or a loop of a 
suitably designed multicore fiber. The FLDW method, 
which has already been applied to a number of 3D quan- 





FIG. 4. Two-photon correlation functions r^'J(nr) (Eq. ( 12 l ) for circular array with no twists (top row) a) n = 0, b) n = 1, 
c) n — 2, d) n = 3 and for the twisted circular array (bottom row) e) n = 0, f) n = 1, g) n = 2, h) n = 3. For clarity we 
re-scale the vertical axis by cos 4(n - 1) (6l)sin 4 (6'). On the twisted array (bottom row) there is a 4 x 27r/12 rotation each loop. 



turn photonic problems [5J E] , can in principle produce 
any of the structures described above including both the 
main loops and the input/output couplers. There are 
some practical challenges in writing waveguides with the 
same properties at different depths within the glass block, 
but these can be significantly ameliorated using wave- 
front correction techniques [17] . On the other hand ex- 
isting multicore fibers [18] could be tapered to induce 
weak coupling between the different cores to create the 
circular and twisted circular arrays. In this approach the 
input/output couplers would be spliced onto the multi- 
core fiber to close the loops. This fiber based approach 
has the additional benefit of large loop radii as fibers of 
lengths of hundreds of metres are readily achievable. 

Since the basic waveguide parameters — index con- 
trast waveguide width, mode size, group index, coupling 
strength, etc, — are similar for these two approaches we 
consider only the FLDW system in detail. We consider 
operation at a wavelength A = 800 nm to exploit stan- 
dard spontaneous down-conversion photon sources and 
efficient silicon avalanche photodiode single photon de- 
tectors, but a configuration for A = 1550 nm is also 
realistic. We next demonstrate that competing limits 
associated with constraints on structure sizes, photon 
loss due to bend-induced radiation, dispersion-induced 
pulse spreading, and coupler bandwidths permit a feasi- 
ble working parameter space. 

Using waveguide parameters comparable to those in 
previous experiments in quantum integrated photon- 
ics [6l 119) , we assume a Gaussian index profile with max- 
imum index contrast Sn — 1.5 x 10~ 3 , a 1/e 2 width 
w = 4.5 /iin and a background index of n^g = 1.44. Typ- 
ical writing setups restrict the total device size to approx- 



imately 40 cm dimensions, so we assume a loop radius of 
R c = 20 cm and length L w 2ttR c . Standard results 
in optical waveguide analysis for Gaussian profile waveg- 
uides [20 then yield a normalized frequency V = 2.32 
which is in the single mode regime when V < 2.405. For 
this loop radius the attenuation coefficient due to bend 
loss is 7 = 6.8 x 10 _7 cm _1 which implies a loss fraction 
of only 0.9% after 100 transits around the loop. 

We suppose the input photons are transform-limited 
with a pulse width of T = 20 ps and a bandwidth 
Auj/(2ir) =50 GHz. The spatial extent of the pulses 
on input is then I = v g T = 0.41 cm. Assuming typi- 
cal silica material dispersion, the net group velocity dis- 
persion is D = —150 ps/nm/km so that the temporal 
width increases by less than 2% after 100 transits around 
the loop. The waveguide dispersion is even weaker and 
may be neglected. Therefore as the photons can propa- 
gate many times around the loop with minimal loss and 
with a spatial pulse length always satisfying l p <C L, it 
is appropriate to treat the photons in the discrete fash- 
ion we have adopted in the analysis above. Furthermore, 
the bandwidth Aw of the pulses corresponds to a wave- 
length bandwidth of AA ~ 17 pm and a relative band- 
width of AA/A = 2.1 x 10~ 5 . Over this bandwidth a 
typical input/output coupler with waveguide separation 
d = 10 /im shows a coupling variation of the same scale 
~2x 10~ 5 , so that it is appropriate to assume a single 
directional coupler reflectivity for the photon pulses. 

Due to the shape of the focal spot of the writing laser 
beam, the waveguides are typically slightly elliptical in 
cross-section. One would arrange to launch the pho- 
tons into one particular polarization axis of the waveg- 
uide, e.g. the vertical. As the waveguides move around 
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each other and the coupling changes from left-right to 
up-down, there will be a small change in the coupling 
strength with position around the loop. However the net 
behavior is still described by a single unitary matrix U 
so there is no significant change to our analysis. 



Finally, we have also assumed the optical path lengths 
in each loop to be the same by tailoring the fan-out sec- 
tion at the directional couplers. Due to fluctuations in 
writing conditions and in the local properties of the glass 
block, this is difficult to accomplish directly. However it 
is possible to post-tune the structures by additional laser 
exposure after characterization and so equalize the path 
lengths. 



Therefore we have shown that the waveguide array de- 
vices considered in this paper are experimentally feasible 
in the near future. In particular we found a parame- 
ter regime where the photon loss due to bend-induced 
radiation, dispersion-induced pulse spreading, and cou- 
pler bandwidths have only a minor effect on the systems' 
dynamics. This justifies the assumptions made in the 
theoretical model in section HH 



VI. CONCLUSION 

In summary, we have studied the dynamical evolution 
of CTQWs on several exotic waveguide arrays in 3D inte- 
grated photonics. Two of the primary benefits of the de- 
sign are the ability to discretely observe CTQWs and the 
ability to inject additional photons during the quantum 
walk. We found that properties of the waveguide array 
topology are readily observable in the time-dependent 
photon correlations. By targeting time steps with par- 
ticular topology induced symmetries or injecting photons 
with a time delay it is possible to engineer a variety of 
multimode entangled states. 
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